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Looking for superconductors with higher transition temperature requires a guiding principle. In 
conventional superconductors, electrons pair up into Cooper pairs via the retarded attraction me¬ 
diated by electron-phonon coupling. Higher-frequency phonon (or smaller atomic mass) leads to 
higher superconducting transition temperature, known as the isotope effect. Furthermore, super¬ 
conductivity is the only instability channel of the metallic normal state. In correlated systems, 
the above simple scenario could be easily violated. The strong local interaction is poorly screened, 
and this conspires with a featured Fermi surface to promote various competing electronic orders, 
such as spin-density-wave, charge-density-wave and unconventional superconductivity. On top of 
the various phases, the effect of electron-phonon coupling is an intriguing issue. Using the func¬ 
tional renormalization group, here we investigated the interplay between the electron correlation and 
electron-phonon coupling in a prototype Hubbard-Holstein model on a square lattice. At half-filling, 
we found spin-density-wave and charge-density-wave phases and the transition between them, while 
no superconducting phase arises. Upon finite doping, d-wave/s-wave superconductivity emerges in 
proximity to spin-density-wave/charge-density-wave phases. Surprisingly, lower-frequency Holstein- 
phonons are either less destructive, or even beneficial, to the various phases, resulting in a negative 
isotope effect. We discuss the underlying mechanism behind and the implications of such anomalous 
effects. 


According to the standard Bardeen-Cooper-Schrieffer 
(BCS) theory of superconductivity, [1] the only possible 
instability of a metallic normal state, described by the 
Landau-Fermi liquid, is the Cooper pairing toward su¬ 
perconductivity (SC) upon an attractive interaction. [2] 
The electron-phonon coupling (EPC) can mediate a re¬ 
tarded attractive interaction between electrons. This has 
to withstand the repulsive Coulomb interaction. Fortu¬ 
nately in conventional metals, the long-range Coulomb 
interaction is well screened and can be effectively re¬ 
placed by a pseudo-potential for quasi-particles below the 
energy scale of the Debye frequency ujd ■ [3, 4] The transi¬ 
tion temperature T c increases (linearly at weak-coupling) 
with uiD, known as the isotope effect. This has been a 
guiding principle in the search of superconductors with 
higher T c , provided that EPC is the pairing glue. How¬ 
ever, the simple BCS scenario could break down in many 
ways. When the Fermi surface touches van Hove singu¬ 
larities, or is nested, density-waves in the spin or charge 
channel would be favorable. A well-known example is 
the Peierls instability toward the charge-density-wave 
(CDW) phase in one-dimensional (ID) electron systems 
with EPC alone. [5] In correlated electrons systems, the 
local interactions are poorly screened, leaving the vari¬ 
ous orders, such as the spin-density-wave (SDW), CDW 
and unconventional SC, close competitors to each other, 
when the Fermi surface is featured with van Hove sin¬ 
gularities and/or nesting. In the Mott limit, the strong 
local interaction leads to the formation of local spin mo¬ 
ments in the first place. The effect of EPC in such cases 
is an intriguing issue. For example, a long standing 
question is whether EPC plays a significant role for d- 
wave pairing in copper-based [6-8] and s±-wave pairing 
in iron-based high-temperature superconductors. [9] As 


a first step toward the issue, one considers theoretically 
a simplest model with local Holstein phonons and lo¬ 
cal Hubbard interactions, the so-called Hubbard-Holstein 
model (HHM). Much effort has been devoted to under¬ 
stand the various orders and the metal-insulator tran¬ 
sition in the HHM in ID [10] and infinite dimensions 
[11, 12]. In view of unconventional SC other than s-wave 
SC (s-SC), such as d-wave SC (d-SC) in layered materials, 
here we consider a HHM on a 2D square lattice. We han¬ 
dle the interplay between electron correlation and EPC 
by the singular-mode functional renormalization group 
(SM-FRG). [13, 14] Our main findings are as follows. At 
half-filling, SDW and CDW competes, but no SC phase 
arises. Upon finite doping, d-SC/s-SC emerges in prox¬ 
imity to SDW/CDW phases. More interestingly, lower- 
frequency Holstein-phonons are either less destructive, 
or even beneficial, to the various phases, resulting in a 
negative isotope effect. 

The 2D HHM is described by the Hamiltonian 
H= -t ^2 ( C L C + h ' C 0 — P X nia ” l " WD X b i bi 

{ij)a iv i 

+U ^(n if - 7))( n 4 - \) +vJ2 nir7 ( b i +bi )’ W 

i icr 

where t is the nearest-neighbor hopping, /r the chemi¬ 
cal potential, U the local Hubbard interaction, lod the 
Holstein phonon frequency, and rj = g /Here 
g is the EPC matrix element and M is the mass of 
the vibrating ion. Henceforth we set t = 1 as the 
unit of energy. The EPC leads to a retarded attrac¬ 
tion n„ = —\WijJ 2 d /{ijj' 2 d + v 2 ), where v is the Mat- 
subara frequency, W = 8 the electron bandwidth, and 
A = g 2 /(MujjjW) an average EPC constant (which de- 
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pends on the spring constant K = Muj 2 d rather than on 
uid independently). 

We treat the correlation effect and EPC by the 
SM-FRG. [13, 14] In a nutshell, the idea is to 
get momentum-resolved pseudo-potential W 234 , as in 
(l/2)c{ 0 .C2 0 ./Vi234C3cr'C4 C r, to act on low-energy fermionic 
degrees of freedom up to a cutoff energy scale A (for 
Matsubara frequency in our case). Henceforth the nu¬ 
merical index labels momentum/position, and we leave 
implicit the momentum conservation/translation symme¬ 
try. Starting from the local U at A = 00 , V can evolve, 
as A is lowered, to be nonlocal and even diverging, due 
to corrections from both U and n, to arbitrary orders 
and in all possible ways (see Supplementary Materials). 
To see the instability channel, we extract from V and n 
the effective interactions in the general CDW/SDW/SC 
channels, 


[ VcDW ](14)(32) = 2[ V + n 0 ] 1234 — [V + Ha ]l243? 

[ VSDW ](13)(42) = — [ V + Ha ] 1234) 

[ VsC ](12)(43) = [ V + Ha ] 1234* (2) 

The left-hand sides are understood as matrices with 
composite indices. Notice that n is local/flat in 
real/momentum space, and n„ = o enters Vcdw because 
n is direct in the charge channel. Since they all orig¬ 
inate from V + n, Vcdw/Sdw/SG have overlaps but 
are naturally treated on equal footing. The divergence 
of the leading attractive (i.e., negative) eigenvalue of 
Vsc/Sdw/Cdw decides the instability channel, the asso¬ 
ciated eigenfunction and collective momentum describe 
the order parameter, and the divergence scale A c is rep¬ 
resentative of T c . More technical details can be found in 
Supplementary Materials. 

Before embarking on full-wedge FRG results, we di¬ 
gress to gain qualitative insights first from an approxi¬ 
mation to FRG. We keep the local part of V only so that 
the FRG reduces to a simple RG. We focus on half filling, 
where the particle-hole symmetry enables us to solve V 
analytically (see Supplementary Materials), 


R + n 0 


(U + n 0 ) exp 


aXW 2 , A ' 

-(1-tan -1 -) 

ujd tt ujd 


,( 3 ) 


where a is a constant of order unity. We find 

(Vsc, Vsdw, Vcdw ) = {V +IIa, — V— n A , V +2n 0 —n A ) 

are bounded, but the behavior of V still provides inter¬ 
esting implications. (1) We observe that Vsc > Vcdw 
for any A > 0, so SC is absent at half filling. In fact, even 
if A = 0 the SC and CDW channels are exactly degener¬ 
ate. (2) If U + n 0 = U — A W = 0, there is a fixed line 
V = 17, on which Vsdw = Vcdw < 0. This implies a 
phase boundary between CDW and SDW. (The local in¬ 
teractions are dispersionless, but the nesting vector (7r, tt) 
decides the CDW/SDW wavevector.) (3) If U — A W > 0 
(or < 0), V is driven more (or less) repulsive so that 
SDW (or CDW) can be enhanced by EPC; (4) A lower 
ui d leads to stronger enhancement of \V + no|, implying 
that softer phonons are beneficial for CDW and SDW 
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FIG. 1. Phase diagrams for (a) ( n) = 1 and (b) (n) = 0.85. 
The color encodes the transition temperature T c . The dashed 
lines enclose the regime in which T c < ICE 4 . Here lod = 0.5. 


in the respective phase regimes, resulting in a negative 
isotope effect for both phases. This effect can also be 
understood from the quasi-particle point of view. By the 
Lang-Firsov approximation (see Supplementary Materi¬ 
als), the bandwidth is narrowed through polaron effect 
by a factor 


z ~ exp 


XW 1 + e^ D 

2 u>d e^ D — 1 


( 4 ) 


where /3 = l/T. On the other hand, U —> U — A W af¬ 
ter the Lang-Firsov transform. The interaction to band¬ 
width ratio becomes r = ( U/W — A )/z. This ratio is 
amplified by 1/z as long as |U — \W\ ^ 0, implying a 
phase boundary U — XW = 0 between CDW and SDW, 
and more interestingly, the amplifying factor 1/z is in 
nice agreement with the exponential factor in Eq. (3), 
provided that A ~ T < wp for the above z to be ap¬ 
plicable. Since we used the bare band as the input, it is 
amazing that our FRG could foresee the effect of pola- 
ronic band narrowing through the renormalization of the 
interaction. The negative isotope effect uncovered above 
is therefore exactly a manifestation of the fact that the 
polaronic effect is stronger for softer phonons. 

We now turn to the full-wedge FRG results. Fig. 1 
is the phase diagram for u>d = 0.5 for two filling levels. 
The color encodes T c versus U and A. The dashed lines 
enclose a regime in which T c < 10~ 4 beyond our interest. 
At half-filling (n) = 1 in Fig. 1(a), there is a phase bound¬ 
ary U = XW separating the CDW and SDW phases. In 
the CDW phase, T c is enhanced with increasing A. In the 
SDW regime, T c exhibits a dome-shaped behavior along 
the A-axis, implying that a weak A also enhances SDW. 























3 



(b) (n)=1, U=3 (d) <n>=0.85, U=3 



FIG. 2. Effects of EPC on the flow of leading eigenvalues 5 
(plotted as 1/5 for clarity) of Vcdw/sdw/sc at A = 1/8 for 
(a) <n> = ! and U = 0, (b) <n> = 1 and U L 3, (c) <„) = 0.85 
and U = 0, and (d) ( n ) = 0.85 and [7 = 3. The phonon 
frequency is indicated in the legend for all panels. For clarity, 
a channel is dropped if it’s |5| is too weak, and 5 is multiplied 
by a factor of 10 in the flow of SC channel in (d). 


These behaviors are exactly what we discussed and un¬ 
derstood in the above simple RG analysis. Moreover, the 
full FRG is able to capture general pairing channels. For 
example, in the SDW phase, Vsc has a negative eigen¬ 
value in the d-wave pairing channel, but it is always less 
diverging than that of Vsdw • This excludes d-SC at half 
filling. The phase-diagram is in full agreement with the 
quantum Monte Carlo (QMC) result [15] on finite-size 
lattices, demonstrating the reliability of our FRG. 

Away from half filling, our result for (n) = 0.85 is 
shown in Fig. 1(b). Since the nesting is no longer at 
the Fermi level, the CDW (SDW) order is stabilized be¬ 
yond a finite threshold A > A c (U > U c ). Below the 
phase boundary of CDW, we find s-SC is established. 
On the other hand, near the SDW phase boundary, d-SC 
emerges. The proximity between these phases is easily 
understood in view of the overlap in the SDW and SC 
channels, and is also known as a manifestation of pairing 
induced by SDW fluctuations. [16, 17] What’s more in¬ 
teresting here is the phase boundaries of both SDW and 
d-SC phases are curvy in the parameter space, imply¬ 
ing that weak (strong) EPC enhances (suppresses) both 
SDW and d-SC. Given the behavior of SDW versus EPC 
we discussed above, however, the anomalous enhance¬ 
ment becomes natural in view of the overlap between 
SDW and d-SC channels. We notice that in an earlier 
FRG work, [18] the EPC (with Holstein-phonons) appears 
to suppress d-SC. We ascribe the difference to the dilute 
frequencies used for II in their case. 

In order to have a closer view of the effects of EPC 
on the various orders, we plot some representative FRG 
flows of the leading eigenvalues S of Vcdw/sdw/sc i n 


Fig. 2 for A = 1/8. Since we are looking for diver¬ 
gence, we drop out IIo and IIa in Eq. (2) to concen¬ 
trate on the flow of the projections of V in the various 
channels. At half-filling with U = 0/3 in Fig. 2(a)/(b), 
the CDW/SDW channel diverges as A is lowered. We 
have checked that for a pure negative-!/ Hubbard model, 
equivalent to uc = oo and U = 0 in (a), the s-SC 
and CDW channels are exactly degenerate, satisfying the 
SO(4) = SU(2) ® SU(2) symmetry, [19] where the excess 
pseudo-SU(2) arises from the particle-hole symmetry at 
half filling. However, a finite wp breaks the pseudo-SU(2) 
symmetry in favor of CDW, [20] since n is a direct inter¬ 
action in the charge channel. For both CDW and SDW 
channels, 1/5 is higher for lower wu, and so is T c . This 
is just the negative isotope effect discussed earlier. 

For (n) = 0.85 in Fig. 2(c)/(d), the CDW/SDW in¬ 
teraction flow at high A is similar to that in (a)/(b) for 
half-filling, since high energy quasi-particles are insensi¬ 
tive to the Fermi level. As A decreases further, however, 
low energy quasiparticles come into play, but the lack of 
nesting limits the phase space for low energy particle-hole 
excitations, so that the SDW/CDW channel eventually 
saturates. In contrary, there is no phase-space restriction 
for Cooper pairing, and upon an attractive pairing inter¬ 
action, either already existing or induced via the overlap 
to CDW/SDW channels, the SC channel is boosted via 
the Cooper mechanism until it diverges. Not surpris¬ 
ingly, we find s-SC/d-SC in relation to the sub-leading 
CDW/SDW channel. More interestingly, the negative 
isotope effect for CDW and SDW clearly also acts on the 
proximiting SC, as is clear in Fig. 2(c) and (d), and this 
is understood as from the channel overlap. The excep¬ 
tion is the case of u>d = 0.1 in (c), which has the lowest 
T c . In fact this is a case in the BCS limit, since A C 1 
and ojd -C W. We shall return to this point below. 

We now check more systematics for the effects of EPC 
on T c of the various phases. For a pure EPC system 
with U = 0, the transition temperature T c is shown in 
Fig. 3(a) and (b). At half-filling in (a), only CDW phase 
is present, and T c clearly drops with increasing o/d, for 
any A. Such a negative isotope effect is discussed above, 
and is in agreement with the result judged from correla¬ 
tion functions measured by QMC on small clusters. [20] 
For (n) = 0.85 in (b), the CDW phase is realized for large 
A, and T c follows the trend in (a) closely. For weaker A, 
the system yields to the s-SC phase. In proximity to 
the CDW phase, we observe that T c for s-SC drops for 
larger uijj. This is understood as caused by the weaken¬ 
ing of CDW fluctuations, so that T c for s-SC eventually 
inherits a negative isotope effect. For even weaker A, 
however, T c for s-SC increases with small u/jj, in a BCS 
fashion. In fact, the FRG reproduces the exact BCS be¬ 
havior T c oc wd for A — / 0 and ujd -C W, if the CDW and 
SDW channels are too weak to affect the pairing chan¬ 
nel (see Supplementary Materials). This is the case for 
A = 1/8 and lod = 0.1 in Fig. 2(c). Therefore, the neg¬ 
ative isotope effect for s-SC occurs only in proximity to 
the CDW phase. 

















(a) (n)=1, U=0 


(a) U=0 


(b) U=3 
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(b) (n)=0.85, U=0 



FIG. 3. T c versus A and ljd with 17 = 0 in (a) and (b), and 
17 = 3 in (c) and (d). The phases are denoted by both text 
and color. The solid lines are equally spaced by AA = 0.025 
in (a) and (b), and A u>d = 0.2 in (c) and (d). 


For a correlated system with (7 = 3, the transition tem¬ 
perature T c is shown in Fig. 3(c) and (d). At half filling 
in (c), a large A drives SDW into CDW, with a phase 
transition at A c = U/W independent of lod- The tran¬ 
sition temperature is always lower for larger u>d, again 
a manifestation of the negative isotope effect. Moreover, 
we observe in (c) a slight enhancement of SDW by a weak 
A and small ojd- This effect is qualitatively explained by 
Eq. (3), and has been discussed previously. To our de¬ 
light, the slight enhancement is consistent with the DCA 
result for U = 8 in Ref. [21]. For the doped case, as 
(n) = 0.85 in Fig. 3(d), the SDW phase yields to d-SC 
phase, and the CDW phase remains for large A. Here T c 
is laterally higher for lower u>d for both CDW and d-SC 
in the respective regimes. A similar case was observed 
but only for the d-SC phase in Ref. [22]. Moreover, in the 
d-SC regime, even though T c decreases with A for ojd > 1, 
it is lifted by a lower uj£> for a given A. Thus lower fre¬ 
quency phonons are at least less destructive to the d-SC. 
On the other hand, there is a marked enhancement of T c 
by a weak A, up to A = 0.2 for ojd = 0.2, which we ascribe 
to the anomalous enhancement of SDW fluctuations as 
revealed in (c). The reason that the negative isotope ef¬ 
fect is observed in the entire d-SC regime is because d-SC 
occurs only in proximity to the SDW phase. 

Finally, we consider the systematics in doping. We 
set ujd = 0.5 for illustration. For U = 0, Fig. 4 (a) 
shows the CDW phase is present at low doping, and s- 
SC at higher doping. For stronger A, a larger doping 
is needed to enter the s-SC phase. In both phases, T c 
decreases with doping, since the density of states p at 
the Fermi level drops. We compare our result to the BCS 




FIG. 4. Doping dependence of T c for c od = 0.5 and (7 = 0 
(a) and (7 = 3 (b). The phases are denoted by both text and 
color. The dashed line in (a) is a fit to the BCS theory (see 
the text for more details). 


formula (dashed line), T^ cs = 1.13u;£>exp(— 1/pVbcs)- 
[1] We choose the value of Vbcs so that T^ cs matches 
our FRG result for A = 0.25 and a deep doping level 1 — 
( n ) = 0.9. We find T c > T^ cs approaching half filling. 
The enhancement follows from the effect of increasing 
CDW fluctuations, also favorable for s-wave pairing but 
missing in the simple BCS theory. For a nonzero (7 = 3 
in Fig. 4(b), the CDW and s-SC phases are realized if A is 
sufficiently large (e.g., A = 0.5). For a weaker A = 0.25, 
d-SC sets in since the SDW fluctuations become stronger. 
Closer to half filling, the SDW phase eventually sets in. 
For both phases, T c is higher for A = 0.25 than that for 
A = 0, reconfirming the previous result that a weak EPC 
enhances SDW/d-SC if wd is small. 

In summary, we investigated the effects of EPC 
in a 2D HHM systematically by SM-FRG. We found 
lower-frequency Holstein-phonons are beneficial to all of 
CDW, SDW, and s-SC/d-SC in proximity to CDW/SDW 
phases, resulting in a negative isotope effect. The quali¬ 
tative mechanism is as follows. For CDW, low frequency 
phonons can be easily softened and adapt to the CDW 
order. Near the phase boundary of CDW, the enhanced 
CDW fluctuations are beneficial to s-SC. For SDW and 
d-SC, the enhancement can be effectively ascribed to po- 
laronic band narrowing, which in turn blows up the corre¬ 
lation effect, favoring SDW and the d-SC in its proximity. 

A few remarks are in order. (1) In a strict 2D sys¬ 
tem, academically there is no finite temperature SDW 
phase by the Mermin-Wagnar theorem, [23] and there is 
only algebraic SC order below the Kosterlitz-Thouless 
temperature. [24] In this regard, T c in our case should be 
understood as a crossover temperature in 2D, or the tran¬ 
sition temperature in quasi-2D systems. (2) We should 
stress that FRG is perturbative in nature, and it works 
best in the itinerant picture up to moderate U/W < 1 
and A < 1 discussed in this paper. We expect applica¬ 
tion to cuprates in the overdoped region, upon necessary 
extension to more realistic phonon modes in cuprates. 
On the other hand, we expect application to the mono- 
layer FeSe on the SrTi03 substrate, where T c is found 
to be much higher than in the bulk. [25] The necessary 
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extension is to account for the multi-orbitals in FeSe. 
Technically the FRG scheme in this paper improves the 
cutoff scheme in Ref. [26], and a further investigation of 
FeSe/SrTi03 is in progress. (3) In the strong correlation 
limit, EPC might also enhance the T c of SDW [21] as a 
result of self-localization of polarons in the presence of a 
sufficiently strong EPC. [21, 27-31]. On the other hand, 
we notice that whether EPC would enhance d-SC in the 
strong correlation limit is under debate. [32, 33] Finally, 
there is even a proposal that EPC-driven bipolarons are 
necessary ingredients for high-T c SC in cuprates. [8] 
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Supplemental Material for “Phonon 
enhancement of electronic orders and 
negative isotope effect in the 
Hubbard-Holstein Model on a square 
lattice” 

SI. TECHNICAL DETAILS OF SM-FRG 

Consider the interaction hamiltonian Hj = 
Vi 234 C 3 (T 'C 4 0 .. Here the numerical index 
labels momentum/position, and we leave implicit the 
momentum conservation/translation symmetry. The 
spin SU(2) symmetry is guaranteed in the above 
convention for Hi. The idea of FRG is to get the 
one-particle-irreducible interaction vertex for fermions 
whose energy/frequency is above a scale A. Equiv¬ 
alently, such an effective interaction is what’s called 
pseudo-potential for fermions whose energy/frequency 
is below A. Starting from the local U at A = oo, the 
contributions to dV/dA are illustrated in Fig. SI. In 
principle there will also be self-energy correction to 
fermions, which we ignore as usual, given the fact that 
we are just looking for the instability of the normal 
state. To proceed, it is useful to define matrix aliases of 
the rank-4 ‘tensor’ V via 

bl234 = -P(12)(43) = C(13)(42) = P(14)(32)- (S5) 

Then dV/dA can be compactly written as 

=[D X ph (V -C) + (T>- C) X ph V](i4 ) (3 2) 

+ [V X PP V} { 12)(43) - [Cx Ph C] {1 3 )( 42), (S6) 

where matrix convolutions are understood within the 
square brackets, and 

p = p + n A , c = c + n A , T) = d + no, 

A P (ab)(cd) = 7 pGa C (A)G bd (-A) + (A -a -A)], 

X P ab)(cd) = ~-^[Gac{A)Gdb{A) + (A -a —A)], (S7) 

where n enters as a matrix (local in real space and flat in 
momentum space), G is the normal state Green’s func¬ 
tion, and we used a hard-cutoff in the continuous Mat- 
subara frequency. Notice that n 0 enters V because the 
EPC induced interaction is direct in the charge chan¬ 
nel. This is also evident from Fig. SI. Since the external 
lines are set at zero frequency (the frequency dependence 
is irrelevant for 4-point interactions in the RG sense), 
the frequency on the phonon lines (thickened wavy lines) 
overlayed by D is automatically zero in Fig. Sl(c)-(e). 

The integration of dV/dA toward decreasing A gener¬ 
ates all one-particle-irreducible corrections to V from U 
and n to arbitrary orders and in all possible ways. We 
extract from V and n the effective interactions in the 
general SC/SDW/CDW channels 

{Vsc, Vsdw,Vcdw) = (P, —C, 2 V - C ). (S8) 
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FIG. SI. One-loop contributions to dV/dA. The greyed bar 
and wavy line denote V and II, respectively. They add up 
where overlayed. Spin is conserved during fermion propa¬ 
gation and is left implicit. The slash denotes the single¬ 
scale propagator and can be put on either one of the fermion 
lines within the loop. The directed-circle indicates circula¬ 
tion of frequency along the loop, and A the running scale. 
The thin (thick) wavy line shares the loop frequency (is at 
zero frequency). Each diagram can be viewed as a convolu¬ 
tion of aliases of V (together with II) via V 1234 = -P(i 2 )( 43 ) = 

C(13)(42) = A>(14)(32). 


(This expression is exactly equivalent to Eq. (2) in the 
main text.) Since they all originate from V + II, they 
are overlapped but are naturally treated on equal foot¬ 
ing. Viewed as scattering amplitude of composite bosons, 
the effective interactions can be decomposed into eigen 
modes. For example, in the SC channel (with a zero 
collective momentum), 

n©d(k,-k)(k',-k') = E fm(k)SM k'), (S9) 

m 

where S m is the eigenvalue, and / m (k) is the eigenfunc¬ 
tion. We look for the most negative eigenvalue, say 
S = min[5©], with an associated eigenfunction /(k). If 
S diverges at a scale A c , it signals the instability of the 
normal state toward a SC state, with a pairing function 
described by /(k). Similar analysis can be performed in 
the CDW/SDW channels, with the only exception that 
in general the collective momentum q in such channels 
is nonzero. Since q is a good quantum number in the 
respective channels, one performs the mode decomposi¬ 
tion at each q. There are multiple modes at each q, but 
we are interested in the globally leading mode among all 
q. In this way one determines both the ordering vec¬ 
tor Q and the structure of the order parameter by the 
leading eigenfunction. Finally, the instability channel is 
determined by comparing the leading eigenvalues in the 
CDW/SDW/SC channels. 

In principle, the above procedure is able to capture 


the most general candidate order parameters. In prac¬ 
tice, however, it is impossible to keep all elements of V 
for computation. Fortunately, the order parameters are 
always local or short-ranged. This is notwithstanding 
the possible long-range correlations between the order 
parameters. For example, the s-wave pairing in the BCS 
theory is local, since the gap function is a constant in 
momentum space. The order parameter in usual Landau 
theories are assumed to be local. The d-wave pairing 
is nonlocal but short-ranged. The usual CDW/SDW or¬ 
ders are ordering of site-local charges/spins. The valence- 
bond order is on-bond but short-ranged. In fact, if the 
order parameter is very nonlocal, it is not likely to be sta¬ 
ble. The idea is, if it is not an instability at the tree level, 
it has to be induced by the overlapping channel. But if 
the induced order parameter is very nonlocal, it must 
be true that the donor channel has already developed 
long-range fluctuations and is ready to order first. These 
considerations suggest that most elements of the ‘tensor’ 
V are irrelevant in the RG sense and can be truncated. 
Eq. (S6) suggests how this can be done. For fermions, all 
4-point interactions are marginal in the RG sense, and 
the only way a marginal operator could become relevant 
is through coherent and repeated scattering in a partic¬ 
ular channel. Therefore, it is sufficient to truncate the 
range between 1 and 2, between 3 and 4, in V{i 2 )(AS)i but 
leaving the range between the two groups arbitrary (thus 
thermodynamical limit is not spoiled). Similar consider¬ 
ations apply to C and T>. Eventually the same type of 
truncations can be applied in the effective interactions 
V CDW/SDW/SC ■ Such truncations keep the potentially 
singular contributions in all channels and their overlaps, 
underlying the key idea of the SM-FRG. [13, 14, 34] The 
merit of SM-FRG is: 1) It guarantees hermiticity of the 
truncated interactions; 2) It is asymptotically exact if the 
truncation range is enlarged; 3) It respects all underly¬ 
ing symmetries, and in particular it respects momentum 
conservation exactly. 4) In systems with multi-orbitals or 
complex unitcell, it is important to keep the momentum 
dependence of the Bloch states, both radial and tangen¬ 
tial to the Fermi surface. This is guaranteed in SM-FRG 
since it works with Green’s functions in the orbital basis. 
We take these as advantages of SM-FRG as compared to 
the patch-FRG applied in the literature. [35-37] 

BCS limit : We notice that if only Fig. Sl(a), the pair¬ 
ing channel, is kept, the BCS theory is trivially repro¬ 
duced. For this to be valid, one requires A c <wd -C W 
and the absence of any nesting, so that the contributions 
from the other channels, Fig. Sl(b)-(e), are negligible. 
To make analytical solution accessible, we approximate 
II„ as a step function, IR = — AW9(ujd — M )- Thus 
IIa = 0 for A > ujd, and the RG flow above 10 d merely 
generates a renormalized Coulomb interaction V*. The 
flow for A < ujd is, with IIa = —AW in the above ap¬ 
proximation, 

d(V- AW)/dA= (p/A)(V- AW) 2 , (S10) 

where p is the normal state density of states, and we 





































assumed that Vk.-k.-k',k' is independent of k and k'. as 
assumed in the BCS theory. (This means that we are 
treating the s-wave pairing channel.) The solution is, 
given the boundary condition at A = cod, 


P-AVF = 


V* - XW 

l + (A-/z*)ln(A /u D y 


(Sll) 


where we used p* = pV* and pW ~ 1. There is a di¬ 
vergence V — A W —> —oo if and only if A — p* > 0 (i.e., 
EPC mediated attraction overwhelms the repulsive V*), 
at the scale 


A c = oj D e~ 1/ix -^\ (S12) 

This is already in nice agreement with the T c in the 
Eliashberg theory, given the approximations in II„. This 
example shows that the idea of pseudopotential can be 
pushed down to any energy scale (not just at ojd as in the 
BCS theory) until it diverges, and the divergence scale 
is just a representative of the transition temperature T c . 
If sufficiently strong, the CDW/SDW channels neglected 
in the BCS theory will clearly invalidate the latter, as 
revealed in the main text. 

Local limit: On the other hand, if only the local ele¬ 
ments of V is kept, we have V = P = C = D. Further¬ 
more, in the presence of particle-hole symmetry (at half 
filling in HHM), the second line of Eq. (S6) cancels out 
(in the local limit), leaving, 

M = ->»-n»)Af(r + n o) . (sis) 

where xa ~ a/A is a local susceptibility at the scale 
A, with a factor a of order unity. This can be solved 


analytically, 


v + n 0 ~ (u + no) exp 

where we used IIo = - A IT. 
text. 


aXW . 2 A v _^ 

-(1-tan X 3 M) 

cod tt wd 


This is Eq. (3) in the main 


S2. POLARONIC BAND NARROWING 


The coupling between electron and phonon can be for¬ 
mally decoupled by the Lang-Firsov transformation. De¬ 
fine a unitary matrix U = exp^/wo) JA rii(bi — fo|)], 
where we recall that rj = g/y/2MuD = \]XWujd /2. It is 
easy to show that H = UHU t becomes 


"H =—t ^2 (cLcjo- + h.c.) - +uj d '22 b l b i 

(ij)a ia i 

+t/]T(n it -l/2)(n 4 -l/2), (S15) 

i 

where Cj = ae~^ ri ^ u ’ D ^ bi ~ b ^ and U = U — XW. The 
hopping part averaged over the phonon ensemble leads 
to a polaronic renormalization of t —> zt, with 


z = exp 


XW 1 + e^ D 
2 u>d e^ D — 1 


(S16) 


This factor describes the coherent part of the kinetic en¬ 
ergy in the presence of EPC, namely the renormalization 
factor for the cohorent bandwidth. For the electrons to 
hop coherently, one requires T -C ojd so that phonon 
excitations are rare, the condition for z to make sense. 













